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Abstract 

' For high energy long baseline neutrino oscillation experiments, we propose a Figure of Merit 

criterion to compare the statistical quality of experiments at various baseline lengths under the 
condition of identical detectors and a given neutrino beam. We take into account all possible 
experimental errors under general consideration. In this way the Figure of Merit is closely 
related to the usual statistical criterion of number of sigmas. We use a realistic neutrino beam 
for an entry level neutrino factory and a possible superbeam from a meson source and a 100 
kt detector for the calculation. We considered in detail four baseline lengths, 300 km, 700 km, 
2100 km and 3000 km, in the neutrino energy range of 0.5-20 GeV for a 20 GeV entry level 
neutrino factory and a 50 GeV superbeam. We found that the very long baselines of 2100 km 
and 3000 km are preferred for the neutrino factory according to the figure of merit criterion. 
Our results also show that, for a neutrino factory, lower primary muon energies such as 20 GeV 
are preferred rather than higher ones such as 30 or 50 GeV. For the superbeam, the combination 
of a long baseline such as 300 km and a very long baseline like 2100 km will form a complete 
measurement of the oscillation parameters besides the CP phase. To measure the CP phase in 
a superbeam, a larger detector (a factor 3 beyond what is considered in this article) and/or a 
higher intensity beam will be needed to put some significant constraints on the size of the CP 
angle. 



1 Introduction 



Neutrino oscillations are to date the only experimental indication of physics beyond the standard 
model that may further expand the horizon of our most basic knowledge of nature. Although 
the existing data from Super-Kamiokande experiment Q and the various other corroborating 
experiments Q offer very strong indications of neutrino oscillations, the oscillation parameters 
have not been determined with sufficient accuracy and an appearance experiment have not been 
convincingly performed. In the ongoing neutrino oscillation experiments and those of the next 
generation under construction, some of the parameters will be probed with greater accuracy 
but others may not yet be accessible. Hence more experiments designed to look for the missing 
information, to carry out appearance experiments, and to probe the known parameters with 
even greater accuracy are desirable. 

In this article we propose a criterion for the determination in general terms and on a global 
scale the overall quality of long baseline (LBL) experiments. We focus on high energy neutrino 
beams in the range of 0.5- 20 GeV with the possibility of r production and CP violation mea- 
surement, and distances from 300 km to 3000 km. At the lower end, actually 295 km baseline, 
is the proposed oscillation experiment ||] using the neutrino beam from the newly approved 
high intensity 50 GeV proton synchrotron in Japan called HIPA || with the Super-Kamiokande 
detector or its updated version. The 700 km baseline is close to the first group of next genera- 
tion LBL experiments: MINOS §, ICARUS @ and OPERA @, all having a baseline of about 
730 km. The 2100 km distance is a possibility for a very long baseline LBL experiment, called 
H2B, under discussion |j, 0], The neutrino beam would be from HIPA and the detector, tenta- 
tively called the Beijing Astrophysics and Neutrino Detector (BAND), will be located in Beijing, 
China. The 3000 km baseline has been discussed extensively in neutrino factory studies |n| 
although particular accelerator and detector sites have not been formally identified. 

To compare different experiments is usually difficult. It requires knowledge of the actual 
detectors, their different systematics, and their neutrino beams. This is beyond the scope of 
the present work. So for a possible yet meaningful comparison we assume that all experiments 
have identical detectors and share a given neutrino beam. We believe that our present idealized 
approach allows us to make a meaningful comparison of the capabilities of the various baseline 
experiments in general terms and will be useful for the choice of the most suitable baseline 
length with a given neutrino beam. 

In Sec. 2 we present some of the fundamentals of LBL experiments which include the two 
types of neutrino beams and their beam fluxes that we will use in our calculations. We also list 
the neutrino charge current cross sections. In Sec. 3 we propose the figure of merit as a criterion 
to determine the quality of an experiment. Section 4 discusses the figure of merit of the various 
measurements in the three-neutrino scheme at four baselines: 300 km, 700 km, 2100 km and 
3000 km. In practice, for the appearance measurement a superbeam will measure — > v ei v r 
or P„ — > D ei D T while a neutrino factory can also measure u e — > Vn , v T or v e — > v r . In order to 
compare the neutrino factory and the neutrino superbeam, in this paper we only use the and 
&a beams in our consideration. The measurements considered include v e and v r appearance, the 
matter effect, leading oscillation parameters, sign of the leading neutrino mass-square difference, 
and the CP violation effect. Section 5 contains a brief summary and some pertinent comments, 
in particular, the comparison of neutrino factories of different primary muon energies. 
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2 Fundamentals of LBL experiments 



Due to the extremely weak interaction cross section of neutrinos and beam divergence at a 
long distance from the source, both the neutrino beam intensity and the detector mass must be 
maximized in order to have the desired statistics. New technologies for both accelerators and 
detectors may be required. In the following we discuss briefly some of the fundamentals of the 
LBL experiment which will be used in our calculations. 



2.1 Neutrino beams 



(2) 



Neutrino factory: The neutrino factory delivers a neutrino beam which contains comparable 
amount of (p^) and v e (z/ e ) obtained from the fi~(p + ) decay in a //-storage ring. 

/!~(/i + ) -> Vptyn) + v e (v e ) + e~(e + ) (1) 

The presence of both muon and electron neutrinos is not a problem because they have opposite 
sign charged leptons in charge current reactions and therefore will not be a background to 
each other if the detector can distinguish between positive and negative electric charges. A 
description of the neutrino factory can be found in JnJ . 

The neutrino flux at baseline L from a neutrino factory of unpolarized muon of energy 
is given by 

d$ _j 2xj(3-2x f )^ ; for u„ 
dE v \ Ux^l-xf)^ for v e 

where Xf = E^/E^, n is the number of useful decaying muons, and 7 = E^/m^ with E^ 
and being respectively the energy and mass of the muon. Two scenarios of the number of 
neutrinos in a beam have been considered: n — 6 x 10 19 /year for an entry level factory and 
no = 6 x 10 20 /year for a high performance factory. For a discussion of the neutrino beam spread 
in a neutrino factor together with some sample plots, see Ref. [O] , 

Meson-neutrino superbeam: Neutrinos from a meson source are obtained from decays 
of pions and kaons produced by collisions of the primary proton with the nuclear target. The 
secondary meson beam produced by the collision is sign selected and then focused by a magnetic 
field. The meson beam is then transported to a vacuum decay pipe, whose length depends on 
the desired energy of the neutrino. The primary neutrino or anti-neutrino beam consists mostly 
of the muon flavor from ^ and K decays in the decay pipe but some impurities of electron 
neutrinos are expected due to a finite branching ratio of , K ± and fi^ decaying into electron 
neutrinos. For example, the NuMI muon neutrino beam at Fermilab contains 0.6% electron 
neutrinos. The neutrino beam energy profile is more complicated. In general, it can be a wide- 
band beam covering a broad range of energies, or narrow band beam with a selected, well-defined 
narrow range of energy. 

The neutrino flux at the detector site is determined by the baseline L, the number of primary 
protons on target (POT), the proton energy E p , and the neutrino energy E v . The following 
empirical formula |13| describes the meson production from a proton beam on a nuclear target: 

da f d 3 a 

XM-, = 27T / XAlEp—rPtdPt 

dx M J dp A 



, 1 4- ^p-Dxm 

2vr / B(l-x M ) A ,J°° TT^j PtdPt 



N M (l-x M ) A (l + 5e- DxM ) (3) 
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where Nm is a normalization factor, p the proton 3-momcntum, xm the Feynman x-variable 
defined as the momentum of the secondary meson divided by the momentum of proton. A, B, 
C, and D are numerical parameters which are different for different secondary particles. Table 1. 
gives their fitted values for 7r + ,7r~,A' + , and AT", taken from Ref.|b|]. By taking the general 
property of Eq. (2), i.e., the E 2 / L 2 behavior, to account for the transverse momentum spread 
due to pion decays, the envelope of the neutrino flux which is the maximum flux at each energy, 
can be written as, 

HE,,,L)« E ^-^l + ^ (4) 

where we take xm — x u = 2E U / E p . The wide or narrow band beam can then be selected from 
this envelope of the total beam. It is interesting to note that the simple formula given in Eq. 
(Q) can account for the beam design of MINOS at various energies as discussed in Ref. PJ. 
We should remark that there is a more complete treatment [|l4| for the energy spectrum of the 
meson-neutrino beam. However the difference with the above expression for E v > 1 GeV is very 
small. Owing to its simpler form, we use the above expression in the following calculations. 





A 


B 


C 


D 


7T+ 


2.4769 


5.6817E-2 


0.57840 


3.0894 




3.5648 


5.0673E-2 


0.68725 


5.0359 


K+ 


1.7573 


6.3674E-3 


0.81771 


5.6915 


K~ 


5.4924 


4.1712E-3 


0.89038 


2.2524 



Table 1: Numerical parameters for meson- neutrino energy spectrum 



As a function of E„ and L the flux of an entry level neutrino factory is given by 

n U) 2 

* f {E v ,L) = 2x 2 (3 - 2z/)^f x 6 x 10 22 (5) 

per thousand tons (kt) per year and per cm 2 , where = 6 x 10 19 /year . We will take the 
primary muon energy to be E v = 20 GeV. 

For the neutrino beam of the meson source we take the flux as 

$ m (E v , L) = (1 - x u ) A {l + 5e- D *»)^nt } X 6 x 10 22 . (6) 

The normalization factor Uq is determined by the total number of v^s in the HIPA broad 
band neutrino beam for E p = 50 GeV at L = 295 km §]. This gives n { Q m) = 13.7 x 10 19 in the 
proper units. We note that Eq. (6) gives only the envelope of the flux, not the actual shape 
of the energy spectrum. A more precise simulation can be made once the exact beam profile is 
given. 

2.2 Interaction cross sections 

The detection of the neutrino flavor is through the charge current interaction. For a neutrino 
energy which is small compared to the mass of the W-boson, the charge current cross sections 
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for the electron and muon neutrino are given by 



"vN — 
a vN — 



0.67 x 10" 38 cm 2 K(GeV) 
0.34 x l{T 38 cm 2 .E„(GeV) 



(7) 
(8) 



For the tau neutrino, the above expression is subject to a threshold suppression. The thresh- 
old for the production of the r is E T = 3.46 GeV. The charge current production cross section 
of the t is usually given numerically as a function of the neutrino energy. We fit the numerical 
cross sections [}L5| from the r threshold to 100 GeV and obtained the following expression 



'uN/"uN 



(E v — Et) 2 
c + CiE v + c 2 El 



6{E V - E T ), 



(9) 



where cq — —84.988, c\ — 18.317, and C2 = 1.194. As shown in Fig. 1, the fit, which is valid 
for E v > 4.0 GeV, is good to within 3%. The difference occurs mostly in the neutrino energy 
region of 20-40 GeV. 




Figure 1: A fit of the numerical ratio of the r to fi charge current production cross 
section from the r threshold to 100 GeV. The boxes are the data points and the dotted 
curve is the fit. The horizontal axis is the neutrino energy in GeV and the vertical axis 
is the ratio of r to \i production cross sections. 



3 Statistics and Figure of Merit 



Let us consider the oscillation: v a — > vp. We define two differential number of particle distribu- 
tions as functions of the neutrino energy and baseline length. In the following we will drop the 
subscript v and use E for the neutrino energy : 



N Qa (E,L) = $(E,L)K t T y a a (E) 
N 0f} (E,L) = <S>(E,L)K t T y a p (E) 



(10) 
(11) 



where K t is in kilotons and T y in years. $(£', L) is the flux of the v a and given in Eqs. (||) 
and (||), <r a (E) and a p{E) are the charge current scattering cross sections of v a and z/g. Then 
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the total number of non-oscillating interacting beam neutrinos at the distance L over an energy 
range is given by 

N 0a (L)= J N aa (E,L)dE, (12) 

The limits of the energy integration depend on the energy range selected which can be over the 
whole energy range of the available beam or a selected bin size. 

The number of charged leptons of flavor (3 appearing in the detector, which is the signal for 
identifying the appearance of up in the v a beam, is given by 

N S p{L) = J P a ^p(E,L)N p(E,L)dE (13) 

where P a —>p(E, L) is the oscillation probability for v a — > vp. Again the integration can be over 
the whole available energy range of the neutrino beam or a bin size. The averaged oscillation 
probability over the given energy range is 

where 

Nop(L) = J N Q p(E,L)dE (15) 

We now define the figure of merit in terms of the physically measurable quantity and its 
uncertainties: 

<P a ^p(L)> 

*M(L) = TJT— (10) 

where S < P a ^p(L) > denotes the uncertainties of < P a ^p(L) >. In the usual definition this is 
the number of sigma that determines the quality of the measurement. There are in general three 
sources of uncertainties: (i) the statistical error in the measurement of the charge lepton of flavor 
(3 which is as usual y/Nsp{L) + fpNop, where fpNop is the number of measured beckground 
events, expressed as the fraction fp of Nop', (ii) the systematic uncertainty in the calculation 
of the number of background events from Nop, which can be denoted as rpfpNop. (iii) the 
systematic uncertainty in the beam flux and the cross section which we denote as gpNsp- The 
total error is the quadrature of all these uncertainties. Then the figure of merit can be written 
as 

F M (L) = 



y/Nsp + fpNop + {gpNspY + (rpfpN p) 2 

< P a ^p(L) > 
\l <Pa lf$L) +h + (9P < P a ^L) >y + {rpfpf 



(17) 



This expression can be generalized to the case when the physical quantity under consideration 
involves two or more averaged probabilities. In that situation, each probability has its associated 
uncertainties. In the uncertainty terms we will refer to the linear terms which are inversely 
proportional to the beam intensity as the statistical terms, and to the quadratic terms as the 
uncertainty terms. 

The value of rp depends on how well the background is calculated. The value of gp depends 
on how well the beam flux and detection cross sections are known. Although fp is not a physical 
quantity, it can be estimated from the background of the beam survival measurement. Some of 
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the error factors for different situations have been estimated in the literature |U5, 17, 1^, [l9L |2C|] . 
Many of the early treatments included only statistical errors; we believe that we have included 
all general uncertainty factors as enumerated above Eq. (O). In Sec. 5 we will compare our 
results with previous analyses. 

The expression of the figure of merit can be simplified if we restrict the energy integration to 
very narrow bins so that the oscillation probability and the charge current (CC) neutrino cross 
sections do not vary significantly within the bin, then P a ^p(E, L) and N p(E, L) are essentially 
constant over the range of integration of the neutrino energy. Hence in this narrow bin (NB) 
approximation we have 

Fm (L)j = = = = g^feS (18) 



where Ej is the central value of the energy of the bin and AE the bin size, both in GeV. 
Generally speaking at shorter baselines when the oscillation probability is small and the number 
of events large, the effect of the background is severe. Typical values of fp for a water Cerenkov 
calorimeter are a few per cent Q . For the value of fp of water Cerenkov ring imaging detector, 
see Ref. g. 

We can define the figure of merit over an energy range for a given baseline length L. It gives 
the overall statistical significance of the measurement over the energy range. Let the energy 
range contains Nf, bins. Then, in the narrow bin approximation, the total figure of merit in the 
energy range is 



]T(Fm(L),) 2 (19) 
\ j=i 



The above expressions can be extended to the case of anti-neutrinos, and that containing 
a mixture of neutrino and anti-neutrino. The quantities associated with the anti-neutrino are 
defined by replacing the neutrino CC cross section by the corresponding anti-neutrino CC cross 
section. 



4 Optimum long baselines 

In this section we apply the Figure of Merit to examine the effectiveness of measurements of 
various relevant quantities at medium and very long baseline oscillation experiments relevant 
to the next generation LBL experiments. We assume identical detectors at all sites with a 
neutrino beam reaching to the different sites. To limit the scope of the analysis, we will only 
discuss the scenario of three neutrinos with MSW-LMA which is the most favorable solution to 
the solar neutrino problem pl[ . The mixing parameters are set as follows: Ato§ 2 ~ Am^i — 
3.3 x 10- 3 eV 2 , Am^i = 5 x lCT^W 2 , tan 2 (9 23 = 1.6, tan 2 6»i 3 = 0.025, and 9 12 = 45°. The 
leptonic CP phase (j> is set to be 90°. The earth density has been chosen to be constant of 
p = 3 g/cm 3 with equal number of protons and neutrons. Then the matter effect constant 
A/E v is 2.3 x lO- 4 eV 2 /GeV. 

As already stated earlier, for the neutrino beam we will take E^ = 20 GeV for the neutrino 
factory and E p — 50 GeV for the meson beam. For the error parameters, we will use rp = 0.1 
and gp = 0.05, independent of the neutrino flavor and beam type. The value of fp, taking the 
value 0.01 and 0.03, may be different for different measurements and neutrino beams and will 
be given explicitly in the various cases considered below. The beam energy considered ranges 
from 0.5 to 20 GeV and the bin size is AE — 0.5 GeV. For very long baselines, the cases of 
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neutrino energy of hundreds of MeV and lower require separate investigations. This is because 
in such lower energy regime, the solar neutrino energy scale will come into play and the Eqs. 
(7) and (8) are on longer valid as the quasi-elastic charge interactions 0] become significant. 
The detector size will be taken to be 100 kt ]23| and the running time 3 years. For simplicity 
we use the narrow bin approximation, Eqs. QlSj ) and ([l9|). We will comment in the Summary 
and Comments section on neutrino factories of different primary muon energies. 

Before we discuss in detail the individual measurements, we summarize the general feature of 
the contributions of the statistical and uncertainty terms. For the neutrino factory the statistical 
terms dominate in the lower energy region whose range increases with the baseline length. Away 
from the dominant region, the statistical and systematic uncertainty terms become comparable 
with the latter being slightly larger. For the meson neutrino beam, the statistical terms are 
dominant in regions at both the lower and upper ends of the beam energies. The two types 
of terms become comparable in the region in between which is sizable. This feature holds 
even when the beam intensity increases or decreases by an order of magnitude from the sample 
meson neutrino beam intensity considered below. Hence the statistical terms are critical in the 
determination of the figure of merit. We can conclude that although the individual figures of 
merit depend on the details of the beam intensity and the various error factors, the relative 
figures of merit of different baseline lengths are not too sensitive to the details of the beam 
profile as long as the correct neutrino beam energy and baseline distance dependence are taken 
into consideration. Hence the relative figures of merit discussed below should be approximately 
valid even when the beam intensity profile is not precisely known. 

For the total figure of merit as a function of baseline, we only consider the neutrino factory. 
The use of Eq. (6) for the meson neutrino beam is not meaningful for the total figure of merit 
in view of the complicated nature of the meson neutrino beam. For comparison we also consider 
the total figure of merit for neutrino factories of 30 and 50 GeV. 

To simplify the numerical calculations, we consider only the case of constant matter density. 
We use the vacuum oscillation formulae for — > is^, v r . For — > v e appearance we use 
the approximate expressions given in Ref. j24j|. For the CP phase sensitivity we use the exact 
formula given in Ref. |2f| . 

4.1 Mixing probability sin 2 26q 3 

The oscillation probability P{v^ — ► v e ) is a direct measurement of sin 2 2#i3 since P(y^ —>■ v e ) °c 
sin 2 2#i3. The statistical significance of P{v^ — > v e ) for a given experiment can be measured by 
the figure of merit. From Eqs. ( |l7j ) and (p"8"|), we have 

F$(L) = <P v ^uAL)> m 

V ^Sr^ + ^„ < P^SL) >) 2 + (r^„J 2 

For the error factor f„ e , contributions from the beam are at the level of 0.6% for meson- 
neutrino beams.^ We conservatively choose f Ve — 0.03 for the conventional meson beam and 
fv e = 0.02 for the neutrino factory. Now the figure of merit can be readily calculated using 
Eqs. ( p0[ ) as functions of the neutrino energy E at four different baselines: 300 km, 700 km, 2100 
km and 3000 km. 

To see the importance of the background let us first consider the figure of merit by ignoring 
the background, i.e., setting r = 0, g = and / = 0. The results are given in Figs. 2(a) and 
(b). As shown in Fig. 2(a), there is an optimum neutrino energy, E opt — 15 GeV, independent 

See Ref. fl] for the water Cerenkov ring imaging detector. 
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of the baseline, where the figure of merit is at a maximum for the neutrino factory. For the 
superbeam, each baseline is associated with a different E opt and the smaller baseline has higher 
figure of merit as shown in Fig. 2(b). In these plots and for all those figure of merit plots that 
follow, the dotted line is for L — 3000 km, the heavy dotted line for L = 2100 km, the dash line 
for L = 700 km, and the dash-dotted line for L — 300 km. For the very long baseline of 2100 
and 3000 km the figures of merit are small and tend to oscillate rapidly at low energies, so we 
do not show them for E < 1.5 GcV. 

We have examined the cases of 30 and 50 GeV neutrino factories and found that their figures 
of merit as a function of the neutrino energy are similar to those of the 20 GcV case. The E opt 
occurs at roughly the same value and the height of the maximum decreases as the primary muon 
energy increases. 

However, the above picture is changed when the effect of experimental uncertainties are 
taken into account. As shown in Figs. 2(c) and (d), the figure of merit is generally reduced. 
For the neutrino factory the figure of merit increases with the baseline length, and for the 
longer distances like 2100 km is much higher than the shorter distance like 300 km. For the 
conventional neutrino beams the reduction of the figure of merit at longer distances is less when 
uncertainties are included, and all distances have comparable figures of merit over the whole 
energy range. Hence the comparison of Figs. 2(c) and (d) with Figs. 2(a) and (b) demonstrates 
the importance of the effect of the background in the determination of the relative merits of 
different baselines. It should be noted that the above results are very sensitive to the value of 
sin 2 2#i3. For a smaller sin 2 2#i 3 , the effect of the background is larger and the longer baselines 
will be better. 

Figure (3) gives the total figures of merit, according to Eq. (19), of three neutrino factories 
as functions of the baseline: 20, 30 and 50 GeV. We extend the baseline to 10 4 km to show the 
variation the total figure of merit. At shorter distance the total figures of merit are roughly 
the same for all three neutrino factories. Their values are small but increase as the baseline 
increases. The factory of lower energy will reach a maximum sooner as the baseline increases. 



For three flavors there is one measurable phase in the neutrino mixing matrix which will give 
rise to a CP- violation effect. Although the effect of CP phase in the hadronic sector is small due 
to the small Jarlskog factor arising from small quark mixing angles, the large or even maximal 
mixing angles of the neutrinos suggest that it is not impossible to have a large effect of leptonic 
CP phase. So we will assume the maximal CP phase, = ir/2. 

The CP phase 4> can be measured by looking at the difference of the oscillation probability 
between Piy^ — ► v e ) and P{v^ — > v e ) or between P(v e — ► j/„) and Piy^ — ► v e ). While the 
former has to have the matter effect removed in order to obtain the CP effect, the later, under 
the assumption of CPT conservation, allows a direct isolation of the matter from the CP effect, 
but it can only be done at a neutrino factory in the case of symmetric matter density which is 
approximately true and identical v e and beam which is generally difficult to arrange. 



The CP phase can be obtained by subtracting the matter effect AP(E, L,<f>) — AP(E, L,0). 
Then the figure of merit for the leptonic CP phase measurement can be written as 



4.2 CP phase 



We define 



AP(E, L, 4>) = P v ^ v . (E, L, 0) - P^ e (E, L, 0) 



(21) 




< AP{<p) > - < AP(0) > 

V( 2 )(0) 
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(a) E(GeV) 






(d) E(GeV) 

Figure 2: Figure of merit for the sin 7 26\^ measurement at a) neutrino factories without back- 
grounds, b) meson-neutrino beams without backgrounds, c) neutrino factories with f=0.02 and 
r=0.1, and d) meson-neutrino beams with f=0.03 and r=0.1. In these figures and those that follow 
the dotted line is for L=3000 km, the heavy dotted line for L=2100 km, the dash line for L=700 
km, and the dash-dotted line for L=300 km. The 2000 and 3000 km curves are not shown below 
1.5 GeV. The vertical axis is the figure of merit given by Eq. ( p0|) and the horizontal axis is the 
neutrino energy. 
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2000 4000 6000 8000 1 -KT 
L (km) 



Figure 3: The total figure of merit for sin 2 2#i3 as a function of baseline length for neutrino factories 
of 20 GeV (dot), 30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline length 
and the vertical axis is the total figure of merit. 



N ' 2 ' w = ( T^^ 1 ^^'^ 1 ' (22 » 

+l£ ^w^ +s2l<PW>2 + <m>2) f- 

where to simplify the notation we have defined all the quantities by retaining the reference to 
the CP phase and dropping the other variables. 

The figures of merit are calculated according to Eq. (22) with / = 0.02 for the neutrino 
factory and / = 0.03 for the conventional beam. Figure [|(a) shows the figure of merit for 
neutrino factories with the full uncertainties for the neutrino factory. For L/E fixed at the 
optimum energy, the sensitivity to the CP phase <f> is a factor of 2 better at L=2100 km than 
that at L=300 km for the neutrino factory. The sensitivity to the CP phase for the conventional 
neutrino beam is shown in Fig. ^(b). A similar optimum energy also exists and the sensitivity 
at L=2100 km is about 20% lower than that at L=300 km at their respective peak values. At 
neutrino factories, the figure of merit over the whole energy range, i.e., the total figure of merit, 
is much larger for 2100 km than for 300 km as shown in Fig. 5. 



4.3 Sign of Am| 2 

Since the vacuum oscillation probability is an even function of Amjy , the presence of the matter 
effect is necessary in order to measure the sign of Am| 2 - The simplest way is to compare 
the measured probability P Ufi ^ Ve with the expected values of P+ ^, Ve and Pj~ ^> Ue , where P + 
assumes Am^ > and P~ for Ato 2 ^ < 0. The channel — > v e may be more advantageous 
to use if Am| 2 < 0. This measurement can only be done if sin 2 2#i3 is sizable and the v e or 
v e appearance signal is statistically significant. Other channels, such as the — > survival 
channel, are less sensitive due to the very small matter dependence. It has been suggested in 
Ref. [ p6| that the muons from all channels at the neutrino factory, e.g. the v e (p e ) — > v^iy^) 
appearance and v^iy^) — > v^(p^) survival channels, be included in a given analysis to reduce 
systematic errors. Here we confine ourselves only to a relatively simple approach to compare 
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(b) E(GeV) 



Figure 4: Figure of merit for CP phase measurement at a) neutrino factories with f=0.02 and r=0.1, 
and b) meson-neutrino beams with f=0.03 and r=0.1. The vertical axis is the figure of merit given 
by Eq. (22) and the horizontal axis is the neutrino energy. 
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Figure 5: Total figure of merit for CP phase measurement for the neutrino factories of 20 GeV 
(dot), 30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline length and the 
vertical axis is the total figure of merit. 
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the relative statistical importance at different baselines. The figure of merit in the present case 
can be written as 

F$(L) = " e " ° (23) 

<P ""7 l%^ +f '" + i9u„ < P^(L) >) 2 + {r v J^Y 



z 




(a) E(GeV) (b) E(GeV) 



Figure 6: Figure of merit for the sign of Am| 2 at a) neutrino factories with f=0.02 and r=0.1, and 
b) meson-neutrino beams with f=0.03 and r=0.1. The vertical axis is the figure of merit given by 
Eq. (|23|) and the horizontal axis is the neutrino energy. 



Figure 6 shows the figure of merit as a function of energy at different baselines for neutrino 
factory and conventional neutrino beams with / = 0.02 for the neutrino factory and / = 0.03 
for the conventional neutrino beam. It is clear from the plots that for both beam types, it is 
much better to have an experiment at longer distance such than at shorter distances. The total 
figure of merit as a function of the baseline length is give in Fig. 7. 



4.4 Matter effect 

The matter effect can be measured by looking at the difference between — > v e and 9^ — > v e , 
although a non- vanishing CP phase, 0^0, can alter the result by as much as 20%. The survival 
channels — > and — > are not very effective for the present purpose since their matter 
effects appear only in nonleading terms. 
The relevant figure of merit is given by 

= <P ^:, <Pl> (24) 



I <Pl>+f I o r 2 f 2 I <P-2> + f I „2( Z- p, >2 I Z- p >2 
V N 0vc (L) + Zr J + No^(L) + 9 l< -n. > +<^2> 



where we have dropped the explicit reference to all variables except where it may cause confusion, 
and denote < Pi >= < P„ -> Ve (L) >, and < P 2 >= < P v ^ Pe (L) >. 
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Figure 7: The total figure of merit for the sign of AttiJ 2 for three neutrino factories of 20 GeV 
(dot), 30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline length and the 
vertical axis is the total figure of merit. 



Figure 8 shows the figure of merit as a function of energy at different baselines for a neutrino 
factory with / = 0.02 and the conventional neutrino beam with / = 0.03. The plots show clearly 
that in both cases it is much better to have an experiment at longer baseline than at a shorter 
one. As show in Fig. 9, the total figures of merit as a function of the baseline length for three 
neutrino factories, 20, 30 and 50 GeV. 

4.5 Precision measurement of Am^ and sin 2 2#23 

Although Am§ 2 anc l sm2 2#23 have been measured at Super-K and hopefully will be further 
improved by K2K, MINOS, OPERA and ICARUS, it is still interesting to measure them at 
different baselines and possibly improve the precision. Am| 2 and sin 2 2# 2 3 can be directly 
related to the survival probability P(y^ — > v^). The figure of merit can be written as 



Figure 10 shows the figure of merit at different baselines for neutrino factories and conven- 
tional neutrino beams. Since the background for muon identification is generally much smaller 
than in the case of the electron, we use / = 0.01 for both the neutrino factory and the con- 
ventional neutrino beam. A distinctive peak corresponding to the value of Am^ 2 L where the 
oscillation is at a maximum can be found and the position of the peak is a good measurement 
of Am| 2 . Here the figure of merit can be large for all energies with longer distances favoring 
measurements at higher energies. The total figure of merit as shown in Fig. 11. 




1- < P > 



(25) 




where < P >=< P v ^ V)i {L) >. 
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(a) E(GeV) 




Figure 8: Figure of merit for the matter effect at a) neutrino factories with f=0.02 and r=0.1, and 
b) meson-neutrino beams with f=0.03 and r=0.1. The vertical axis is the figure of merit given by 
Eq. ((23) and the horizontal axis is the neutrino energy. 
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Figure 9: The total figure of merit of the matter effect for three neutrino factories: 20 GeV (dot), 
30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline length and the vertical 
axis is the total figure of merit. 
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Figure 10: Figure of merit for Ami and sin 2 2#23 at a) neutrino factories with f=0.01 
and b) meson-neutrino beams with f=0.01and r=0.1. The vertical axis is the figure of 
by Eq. ([25]) and the horizontal axis is the neutrino energy. 
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Figure 11: The total figure of merit for Am^ and sin 2 2#23 for the three neutrino factories: 20 GeV 
(dot), 30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline length and the 
vertical axis is the total figure of merit. 
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4.6 — >• i> r Appearance 

The appearance of v r from a i/„ beam is an unambiguous signal of the v„ — ► ^ r oscillation. 
Precision measurement of this oscillation probability is crucial in establishing the oscillation 
patterns and in determining whether or not ~ > 2/ T is the only dominant oscillation mode or if 
there is still room for the — > ^ s oscillation, where ^ s is a sterile neutrino. Although we expect 
that in the next few years, K2K and/or OPERA will observe the production of r because of 
the large — > v T probability, it is desirable that future detectors possess good r identification 
capability. The figure of merit can be written as 

4f(D = t= , r ^< P — AL)> (26) 



\I <P ^AL) >+S + r ^ + (g < ^.^W » 2 



Figure 12 shows this figure of merit as a function of energy at different baselines. It is 
expected that the error in the present case will be larger than measurements of other physical 
variables due to r event selection. Hence we set f=0.03 for both neutrino factories and meson- 
neutrino beams. Due to the threshold of tau production, we made a cutoff at E v = 4 GeV. For 
the neutrino factory the larger distance is always better. For the meson beam 2100 km is also 
better than 300 km for neutrino energy above 6 GeV. 

The total figure of merit is shown in Fig. 13. As shown, for L > 1000 km the total figure 
of merit favors factory of higher energy. However, since the figure of merit of all factories are 
high, it is not necessary to maximize the figure of merit. A 20 GeV neutrino factory is more 
than adequate. 





12 16 20 
(b) E (GeV) 



Figure 12: Figure of merit for tau appearance at a) neutrino factories with f=0.03 and r=0.1, i 
b) meson-neutrino beams with f=0.03 and r=0.1. The vertical axis is the figure of merit given 
Eq. (pS) and the horizontal axis is the neutrino energy. 



17 



32 



i 



J I I L 



2000 



4000 6000 
L (km) 



8000 1 -10 



Figure 13: The total figure of merit for the r appearance measurement for the three neutrino 
factories: 20 GeV (dot), 30 GeV (dash), and 50 GeV (dashdot). The horizontal axis is the baseline 
length and the vertical axis is the total figure of merit. 



5 Summary and Comments 

In Table 2 we list the maximal values of all the figures of merit for all the measurements 
considered in the range 0.5 to 20 GeV for a neutrino factory of 20 GeV muon and a superbeam 
of 50 GeV primary proton. To put things in perspective, a few remarks are due. The value of 
the figure of merit is equal to the same value of a in statistical significance. For a meaningful 
measurement we need to have a figure of merit no less than 3. Since the values of the figure of 
merit are obtained in this paper with explicit assumptions about the neutrino flux, the detector 
size, and the running time, we cannot take the values in the table too literally. However, as 
discussed in the preceding section, the relative merits of different baseline, which are not too 
sensitive to the explicit experimental conditions, are meaningful. The table shows that a longer 
baseline is preferred over a shorter one for the neutrino factory. For the meson neutrino beams, 
the preference depends on the variables measured but longer baseline length is generally slightly 
more preferred. For a more complete measurement of all the relevant variables, the combination 
of a relatively short and long baselines, such as 300 and 2100 km, will do a much better job. 

The CP phase measurement is a challenging case. The values of the associated figure of 
merit at individual energies are generally small in all the cases we have considered. This will 
be true even for the total figure of merit in the case of superbeam. For the neutrino factory the 
value of the total figure of merit is already reasonably large at 2100 km for the 20 GeV factory. 
Generally speaking it is desirable to improve the statistics by improving the running conditions 
in all aspects: the beam luminosity, the running time, and the detector size. 

Since the flux of the neutrino factory given in Eq. (5) is realistic, it is appropriate to make 
some detailed observation in the case of neutrino factory on the dependence of the figure of 
merit on the primary muon energy and baseline length: 

(a) First consider the figure of merit at a fixed distance as a function of the neutrino energy. 

From the neutrino flux of the neutrino factory, Eq. (pf), it can be seen that for a given neutrino 

( f) 

energy and n , the neutrino flux decreases with increasing for E^ > AE v /3. Therefore at a 
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neutrino factory 


meson-neutrino beam 




300 km 


2100 km 


300 km 


2100 km 


sin'' 26 13 


6.0 


13 


5.2 


3.9 


CP phase S 


0.55 


1.1 


0.6 


0.5 


sign of Am§ 2 


0.7 


8.0 


1.0 


3.2 


matter effects 


0.5 


6.6 


0.6 


2.4 


Am§ 2 and sin 2 26*23 


75 


165 


95 


58 


t appearance 


9 


17.5 


4.2 


4.7 



Table 2: Summary of relative figures of merit (at the maxima) for various measurements 
at the baselines L=300 and 2100 km. 



given baseline, the maxima of the figures of merit of the various measurements moves to higher 
neutrino energies, in all cases very slowly, when E^ increases, while the height of the maximum 
generally decreases significantly. This is born out in our investigation of the 30 GeV and 50 
GeV factories. 

(b) Next consider the total figure of merit as a function of the baseline obtained by integrating 
the figure of merit over the whole available energy range at a fixed baseline. The total figure of 
merit increases almost linearly with baseline, starting at small baseline, with the values of the 
total figure of merit about the same for all three factories. The values of the figure of merit of the 
three factories will separate with further increases in the baseline. The lower energy factory will 
reach a maximum first. All the maxima are very broad. The first maximum in all the variables 
is reached by the 20 GeV factory at about 2000 km for the tau appearance experiment. 

(c) All the values of the total figure of merit for the neutrino factories we considered are very 
significant, larger than 10, for baseline beyond 2000 km, except for the CP phase. The figure of 
merit for the CP phase of 90° has maxima in the 3500 km and 4000 km region. The height of 
the maximum decreases with increasing neutrino factory energy. The 20 GeV factory has the 
maximum figure of merit value of 4.5 at 3500 km and at 2000 km it is 3.6, which are reasonably 
large. 

(d) Based on the above two considerations the lower energy factory of 20 GeV is preferable 
considering its lower budget of construction and possibly higher physics output. However for 
a more complete investigation, a detailed simulation on the actual number of events should be 
performed. 

(e) All of the studies on neutrino factories use similar principles in determining the best 
L and E combination, i.e., there is some goodness of fit criterion. The differences are in the 
details, such as which sources of uncertainties are included, the particular values assigned, and 
whether or not error correlations are taken into account. Our total figure of merit for neutrino 
factories gives results similar to those obtained in other studies. For example, we find that the 
optimum L for detecting a nonzero sin 2 20\ 3 in a neutrino factory is about 6000-7000 km, and 
that higher stored muon energies do better; similar conclusions are reached in Ref. |l^, pb] |. 
The optimum L for detection of CP violation is 3000-4000 km, which agrees with similar 
studies jl6|, [IT], [l8|, [l9| p0| . We find that the optimum L for determining the sign of Am| 2 is 



around 5000 km, consistent with the conclusions of Refs. 17 and [|2(J. However, we have also 



systematically applied the figure of merit to other oscillation variables, such as the matter effect, 
v T appearance and 023, most of which have not been applied before. Furthermore, we have made 
a comparison of neutrino factories and meson neutrino beams using identical criteria. 

Although we have not considered the total figure of merit for the case of the meson neutrino 
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due to the lack of generally valid beam energy profile, we can nevertheless draw some useful 
conclusions based on the salient features of the differential figure of merit. The maximum of the 
figure of merit at a shorter baseline tend to be sharply peaked, while that of a longer baseline 
is broad. Hence the total figure of merit integrated around the maximum will favor longer 
baselines. Naively integrating the differential figure of merit over the envelope of the beam 
intensity, Eq. (6), shows that the values at 2100 km are larger than at 300 km. 

Finally, we note that in the present work, we focus on the statistical significance of long 
and very long baselines in the energy range of 0.5-20 GeV. Consequently the dominant mass 
scale is the atmospheric neutrino scale, i.e., Am§ 2 - For neutrino energies of a few hundred 
MeV or smaller, the effect of the solar neutrino scale becomes important. Therefore a separate 
investigation is necessary. 
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